Effects of a phase transition on HBT correlations in an integrated 
Boltzmann+Hydrodynamics approach 
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A systematic study of HBT radii of pions, produced in heavy ion collisions in the intermedi- 
ate energy regime (SPS), from an integrated (3+l)d Boltzmann+hydrodynamics approach is pre- 
sented. The calculations in this hybrid approach, incorporating an hydrodynamic stage into the 
Ultra-relativistic Quantum Molecular Dynamics transport model, allow for a comparison of differ- 
ent equations of state retaining the same initial conditions and final freeze-out. The results are also 
compared to the pure cascade transport model calculations in the context of the available data. 
Furthermore, the effect of different treatments of the hydrodynamic freeze-out procedure on the 
HBT radii are investigated. It is found that the HBT radii are essentially insensitive to the details 
of the freeze-out prescription as long as the final hadronic interactions in the cascade are taken into 
account. The HBT radii Rl and Ro and the Ro/Rs ratio are sensitive to the EoS that is employed 
during the hydrodynamic evolution. We conclude that the increased lifetime in case of a phase 
transition to a QGP (via a Bag Model equation of state) is not supported by the available data. 
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I. INTRODUCTION 

One of the main purposes of the research in heavy ion 
collisions (HICs) at high beam energies is to explore the 
existence of the quark gluon plasma (QGP) as well as its 
properties. The equation of state (EoS) of nuclear matter 
is one of the key points to gain further understanding 
since the EoS directly provides the relationship between 
the pressure and the energy at a given net-baryon density. 
Phase transitions (PT), e.g., from the hadron resonance 
gas phase (HG) to the color-deconfined QGP (see e.g., [l|, 
Si El) j constitute themselves in changes of the underlying 
EoS. 

Although, on the low temperature side (and for low 
baryo-chemical potentials /is), investigations of the EoS 
of nuclear matter have been pursued for many years 
and uncertainties have been largely reduced, on the high 
temperature side, the EoS of hot and dense QCD mat- 
ter is still not precisely known. For systems created 
in the RHIC energy region with high temperatures and 
low baryo-chemical potential, lattice quantum chromo- 
dynamics (1QCD) (see, e.g., Ref. Q) calculations predict 
a cross-over transition between the hadron gas and the 
QGP phase. The additional structures of the phase di- 
agram are still under heavy debate, especially regarding 
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the existence or non-existence of a critical endpoint [5[. 

The intermediate SPS energy regime still raises a lot of 
interest because the onset of deconfinement is expected 
to occur at those energies and the possibility of a critical 
endpoint and a first-order phase transition is not yet ex- 
cluded. Several beam-energy dependent observables such 
as the particle ratios!, Q, the flow d & G3, the HBT 
parameters 

[HElEah show 

a non-monotonic behaviour 
around E\ a b = 30 — 40A GeV and the interpretation re- 
mains still unclear. Therefore, future energy scan pro- 
grams at RHIC, SPS and FAIR are planned to explore 
the high-/is region of the phase diagram in more detail. 

To learn something about the hot and dense stage of 
the collision from the final state particle distributions, 
a dynamical modeling of the whole process is necessary. 
Some of the important ingredients which have to be con- 
sidered in a consistent manner are 

• the initial conditions and the initial nonequilibrium 
dynamics, 

• the treatment of the phase transition and 
hadronization, as well as the right degrees of free- 
dom, 

• viscosity effects in the initial partonic as well as in 
the hadronic stage of the evolution, 

• hadronic rescatterings and freeze-out dynamics. 

We notice that part of these have been pointed out to 
be of importance esp ecially for the understanding of the 
HBT results IS Eft- 
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Combined microscopic+macroscopic approaches are 
among the most successful ideas for the modeling of the 
bulk properties of HICs 0, E, d. Recently, a trans- 
port approach that embeds a full (3+1) dimensional ideal 
relativistic one fluid evolution for the hot and dense stage 
of the reaction has been developed and first results are 
convincing [2^, S3 . This hybrid model inherits the ad- 
vantages of the Ultra-relativistic Quantum Molecular Dy- 
namics (UrQMD) model for the dynamic treatment of the 
initial and the final state by taking into account event- 
by-event fluctuations. Furthermore, the hybrid model 
allows for a dynamical coupling between hydrodynam- 
ics and transport calculation in such a way that one can 
compare calculations with various EoS during the hydro- 
dynamic evolution and with the pure cascade calculations 
within the same framework. 

It is well-known that by using HBT interferometry 
techniques one can get detailed information about the 
space-time configuration of the system at freeze-out. We 
concentrate here on the two (identical) pion interferom- 
etry and test the sensitivity of the HBT results on dif- 
ferent stages of the evolution. In our previous investi- 
gations on the HBT correlation of various identical par- 
ticle pairs from HICs at AGS, SPS, and RHIC energies 
[IlSlSiSISESI] , we adopted the UrQMD model but 
further considered the mean field potentials for both con- 
fined and "pre-formed" particles in the model [l3|,SES|]- 
We found that also initial stage interaction may con- 
tribute to a better description of the HBT time-related 
puzzle throughout the energies from AGS, SPS, up to 
RHIC. 

In this paper we perform a systematic investigation of 
the sensitivity of HBT correlation of negatively charged 
pions to the EoS by applying the newly developed hybrid 
approach. Similar more focused studies were frequently 
discussed with microscopic transport or hydrodynamic 
models before [ll|, S3, SM SH- R is a l so interesting to 
study if the current set of EoS employed in the hydro- 
dynamic phase support the conclusion about the origin 
of the HBT time-related puzzle. In addition, the ef- 
fects of the hadronic rescattering and of resonance decays 
(dubbed as "HR") after the hydrodynamic freeze-out on 
the HBT radii and the Ro / Rs ratio deserve more investi- 
gation. We have also noticed that some recent progresses 
of this topic both from an improved hydrodynamic cal- 
culation [16[ and from the pion-optical-potential point of 
view [30| have been published which provides additional 
new insights. 

The paper is arranged as follows. In the next section, 
the UrQMD+hydrodynamics hybrid model is introduced 
briefly. The set of different EoS that are employed in 
the hydrodynamic phase are explained. Two different 
treatments for the transition process from the hydrody- 
namic evolution to the final state hadronic cascade are 
discussed. In Section 3, the analyzing program CRAB 
for constructing the HBT correlator and the correspond- 
ing three dimensional (3D) Gaussian fitting process are 
introduced. In Section 4, the HBT radii Rl, Rq, and 



Rs, and the Ro/ Rs ratio of the negatively charged pion 
source from central Pb+Pb collisions at SPS energies are 
shown and discussed in the context of the experimental 
data. Finally, in Section 5, a summary and an outlook 
are given. 



II. URQMD +H YD RO DYNAMIC MODEL 

An integrated Boltzmann+hydrodynamics transport 
approach is applied to simulate the dynamics of the heavy 
ion collision. To mimic experimental conditions as realis- 
tic as possible the initial conditions and the final hadronic 
freeze-out are calculated using the UrQMD approach. 
Especially for an observable like HBT radii it is impor- 
tant to take care of the complexity of the different effects 
[T3.[l5|. The non-equilibrium dynamics, e.g. fluctuations 
of the local baryon and energy density [3l| , in the very 
early stage of the collision and the final state hadronic 
interactions are properly taken into account on an event- 
by-event-basis. 

UrQMD is a microscopic transport approach based on 
the covariant propagation of constituent (anti-)quarks 
and diquarks accompanied by mesons and baryons, as 
well as the corresponding anti-particles, i.e., full baryon- 
antibaryon symmetry is included. It simulates multi- 
ple interactions of ingoing and newly produced parti- 
cles, the excitation and fragmentation of color strings 
[32L l33l 1341 and the formation and decay of hadronic res- 
onances [35, 36] . In principle it is also possible to incorpo- 
rate mean field interactions in the transport calculation, 
In the present calculation they are neglected in order to 
test the hydro-phase and in the following we will refer to 
the pure cascade calculation as UrQMD-2.3. Studies on 
the thermodynamic properties of UrQMD can be found 
in [iimili • 

The coupling between the UrQMD initial state and 
the hydrodynamical evolution proceeds when the two 
Lorentz-contracted nuclei have passed through each 
other, t start = 2i?/^ 7 2 - 1 [20]. After that, a full (3+1) 
dimensional ideal hydrodynamic evolution is performed 
using the SHASTA algorithm [J, [Io[ . The hydrodynamic 
evolution is stopped, if the energy density e of all cells 
drops below five (default value) times the ground state 
energy density £0 (i.e. ~ 730MeV/fm ). This criterion 
corresponds to a T-/^-configuration where the phase 
transition is expected - approximately T = 170 MeV at 
fiB = 0. The hydrodynamic fields are then mapped to 
particle degrees of freedom via the Cooper- Frye equation 
on an isochronous (in the computational frame) hyper- 
surface. The particle vector information is then trans- 
ferred back to the UrQMD model, where rescatterings 
and final decays are calculated using the hadronic cas- 
cade. We will further refer to this kind of freeze-out pro- 
cedure as the isochronous freeze-out (IF) . This procedure 
is explained in detail in 21 1. 

In this paper we introduce another freeze-out proce- 
dure to account for the large time dilatation that occurs 
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for fluid elements at large rapidities. Faster fluid elements 
need a longer time to cool down to the same temperatures 
than the cells at midrapidity since the hydrodynamic cal- 
culation is performed in the center-of-mass frame of the 
collision. At higher energies the isochronous hypersurface 
increasingly differs from an iso-r hypersurface (r is the 
proper time). To mimic an iso-r hypersurface we there- 
fore freeze out transverse slices, of thickness Az = 0.2fm, 
whenever all cells of that slice fulfill our freeze-out crite- 
rion. For each slice we apply the isochronous procedure 
described above separately. By doing this we obtain a 
rapidity independent freeze-out temperature even for the 
highest beam energies. For lower energies (i?; a fc ^ &0A 
GeV) the two procedures yield very similar results for 
the temperature distributions. The hydrodynamic fields 
are then again mapped to particle degrees of freedom via 
the Cooper-Frye equation on this new hypersurface. In 
the following we will refer to this procedure as "grad- 
ual freeze-out" (GF). A more detailed description of the 
hybrid model including parameter tests and results for 
multiplicities and spectra can be found in 21 1 . 

Serving as an input for the hydrodynamical calcula- 
tion the EoS strongly influences the dynamics of an ex- 
panding system. In this work we use three different EoS 
to investigate their effect on the extracted HBT radii. 
The first EoS, named the hadron gas (HG), describes a 
non-interacting gas of free hadrons [4l|. Included here 
are all reliably known hadrons with masses up to 2 GeV, 
which is equivalent to the active degrees of freedom of the 
UrQMD model (note that this EoS does not contain any 
form of phase transition) . This purely hadronic calcula- 
tion serves as a baseline calculation to explore the effects 
of the change in the underlying dynamics - pure trans- 
port vs. hydrodynamic calculation. The second EoS, 
named the Bag Model EoS (BM), follows from coupling 
a bag model of massless quarks and gluons to a Walecka 
type of hadron gas including only SU(2) flavours (for de- 
tails the reader is referred to [l(). This EoS exhibits 
a strong first-order phase transition (with large latent 
heat) for all baryonic chemical potentials /jb- The third 
EoS, named the chiral+HG (CH), follows from a chiral 
hadronic SU(3) Lagrangian and incorporates the com- 
plete set of baryons from the lowest flavour-S'C/(3) octet, 
as well as the entire multiplets of scalar, pseudo-scalar, 
vector and axial- vector mesons [42[ . Additional baryonic 
degrees of freedom are included to produce a first-order 
phase transition in certain regimes of the T-fj, q plane, de- 
pending on the couplings [H, [43|, |44[. Using this EoS, 
a phase structure including a first-order phase transition 
and a critical endpoint at finite /is is obtained [45]. This 
EoS has already been successfully applied to a hydrody- 
namic calculation [20| . 

To visualize the differences of these EoS, Fig. Q] shows 
the average pressure of the expanding system, from cen- 
tral Pb+Pb collisions at E lab = 20A GeV (left plot) and 
158A GeV (right plot), as a function of time (in the cen- 
ter of mass frame). The vertical line in each plot indi- 
cates the starting time of the hydro evolution. The mean 
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FIG. 1: Time (in the center of mass frame) evolution of the 
average pressure for all three EoS and a central Pb+Pb col- 
lision at _Ei ab = 20A GeV (left plot) and Si ab = 158^4 GeV 
(right plot). The vertical line in each plot indicates the start- 
ing time of the hydro evolution. 



value of the pressure has been obtained by weighting the 
pressure, Pj, in every cell i by its energy density, Ej, and 
integrating over the hydrodynamic grid 
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All curves in Fig. [T] arc plotted until the point in 
time when the isochronous freeze-out criterion is fulfilled. 
Compared to the HG the BM-EoS leads to a delayed 
freeze-out time (i.e. a much longer expansion). While 
in the first few fm/c of the evolution, the system obey- 
ing the BM-EoS expands most violently (due to the high 
pressure gradient in the QGP phase) , once the system en- 
ters the mixed phase, its expansion is slowed down con- 
siderably. This can be observed as the "kink" in Fig. 
[U At the higher beam energy, Ei a b = 158 A GeV, this 
softening of the EoS is even more pronounced. Since 
the HG-EoS does not contain any phase transition, no 
softening can be observed, resulting in the shortest ex- 
pansion time. The chiral CH-EoS lies in between both 
extreme cases. Although a small kink can be observed, 
it is not as pronounced as in the BM-EoS. The effect of 
changes in the EoS on HBT results has been studied be- 
fore [28|], but the great advantages of our approach are 
the full (3+1) dimensions and the same initial conditions 
and freeze-out for all three cases and all beam energies 
without adjusting additional parameters. 



III. CRAB ANALYZING PROGRAM AND THE 
FITTING PROCESS 

To calculate the two- particle correlator, the CRAB 
program is adopted [461 ] , which is based on the formula: 

£,/ k q ) = I d 4 x i d 4 x J g(x l ,p i )g(x j ,p J )\(j)(r',q')\ 2 
J d 4 Xig(xi,pi)J d i x j g(xj,p j ) 

Here g(xi,pi) is an effective probability for emitting a 
particle i with 4- momentum pi = p i ) from the space- 
time point Xi = (ti,Ti). 0(r',q') is the relative wave 
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function with r' being the relative position in the pair's 
rest frame, q = p; — Pj and k = (p i + Pj)/2 are the 
relative momentum and the average momentum of the 
two particles i and j. 

In this work, we select central (< 7.2% of the to- 
tal cross section <jt) Pb+Pb collisions at SPS energies: 
E b = 20A, 30A, AOA, 80A and 158A GeV, with a pair 
rapidity cut < 0.5 (Y^ = log((£i + E 2 + p\\i + 

P\\2)/{Ei + E2 — p\\i — p||2))/2 is the pair rapidity with 
pion energies E\ and E2 and longitudinal momenta p\\i 
and p|| 2 in the center of mass system). For each EoS 
about 2500 events are calculated. All particles with 
their phase space coordinates at freeze-out are then given 
into the CRAB analyzing program. Only the negatively 
charged pions are considered during the analyzing pro- 
cess (for each analysis, one hundred million pion pairs 
are considered). For the cascade calculations, we take 
the results from our previous publications as reference 
[l3l [23j . We found that the residual Coulomb effect af- 
ter the hadron freeze-out on the HBT radii of the pion 
source is small 25], therefore we omit it in the present 
analysis. Finally, we choose the longitudinal comoving 
system (LCMS) frame of the pair (also called the "Out- 
Side-Long" system, in which the longitudinal component 
of the pair velocity vanishes), which is frequently adopted 
in recent years, and fit the correlator by a 3D Gaussian 
distribution 



C(q ,qs,qL) =K[l + \ 

xeM-Rlti ~ RWo ~ RWs ~ 2R 2 OL q oqL )\. (3) 

Here K is the overall normalization factor, the q x and 
R x are the components of the pair relative momentum 
and homogeneity length (HBT radius) in the x direction, 
respectively. The A parameter is called as the incoherence 
factor or, more correctly, the intercept parameter and lies 
for Bose-Einstein statistics between and 1 for two-boson 
correlations in realistic HICs. Because the parameter A 
might be influenced by many additional factors, such as 
contamination, long-lived resonances, or the details of 
the residual Coulomb modification, we regard it as a free 
parameter and do not show it in this letter. However, 
In [25| we have found that the calculated A factor with 
UrQMD can be comparable with experimental data at 
RHIC energies (although somewhat larger than data). 
At SPS energies, ArjrQMD ~ 0.8 — 0.9, is also compatible 
with experimental data which Adata ~ 0.6 — 0.8. The 
R OL represents the cross-term and plays a role at large 
rapidity. To fit the correlator with Eq. (J3J), we use ROOT 
|47j software and minimize \ 2 - 



IV. HBT RESULTS 

Fig. [2] shows the transverse momentum kx (k-x — 
(Pit+P2t)/2) dependence of the HBT radii Rx, Ro, and 
Rs (at midrapidity) of 7r~ source from central Pb+Pb 
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FIG. 2: Transverse momentum kx dependence of the HBT 
radii Rl, Ro, and -Rs (at midrapidity) of n~ source from 
central HICs at SPS energies (£ lab = 20^4, 30A, 40A, 80A, 
and 158A GeV). The NA49 data are indicated by solid stars 
[4^ |. The pure cascade calculation is depicted by lines while 
the hybrid model calculations with different EoS (HG, BM 
and CH) are depicted by dashed lines with open symbols. 
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FIG. 3: Freeze-out time dependence of the it emission in 
central Pb+Pb at E lab = 20.4 GeV (left plot) and 158A GeV 
(right). Calculations with the UrQMD cascade are compared 
with the hybrid model calculations with the EoS of HG, BM, 
and CH. 



collisions at SPS energies. The data (solid stars) are from 
the NA49 Collaboration [48T |. The pure cascade calcula- 
tion is depicted by lines while the hybrid model calcula- 
tions with different EoS (HG, BM and CH) are depicted 
by dashed lines with open symbols. As was shown be- 
fore, the cascade calculation gives a fairly good result of 
the kr dependence of Rl and Rs values except at quite 
small kx, while for Ro, it is slightly larger than data 
at large kx- In contrast, the hybrid model calculations 
show large HBT (for all employed EoS, but to a vary- 
ing degree) in all directions, especially in the longitudi- 
nal direction. The HG and CH are moderately increased 
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and lead to very similar results for all three directions. 
The large latent heat in the bag model leads to a further 
strong increase in the longitudinal direction and in the 
transverse direction at large fey. This increase in the BM 
mode becomes more pronounced at higher beam energies. 
At first glance, this result might be surprising because at 
least in the transverse direction one would expect a faster 
expansion including a hydrodynamic evolution. On the 
other hand, one knows that the system spends a longer 
time without emitting any particles in the hybrid model 
calculation. 

Fig. [3] exhibits the freeze-out time dependence of the 
7T~ emission in central Pb+Pb at E\ a \> = 20A GeV (left 
plot) and i?i a b = 158A GeV (right plot). It is clearly seen 
that there are almost no pions emitted before ~ lOfm/c 
in the hybrid model calculations. This is easy to under- 
stand because even in the gradual hydro- freeze-out which 
is applied here, it takes a while until the first slices have 
cooled down and are frozen out from the hydrodynamic 
evolution. There is no particle emission from earlier times 
in contrast to the pure cascade calculation. For the BM- 
EoS, this effect is present even for a longer time since the 
expansion lasts longer [52]. The on-the-average longer 
freeze-out time leads to an apparently larger size of the 
pion source. Furthermore, it is clear (and expected) that 
the EoS with larger latent heat (such as in the BM mode) 
leads to a longer emission duration of the particles (as 
seen in Fig. [3] when tf > 15fm/e) so that it produces 
larger HBT radii. This behaviour is clearly seen in the 
more time-dependent directions Rl and Ro- This be- 
haviour might be improved by allowing particles also to 
freeze out and fly into the detector at all times of the col- 
lision. Especially, if there are fast pions produced during 
the first hard collisions in UrQMD at the edge of the sys- 
tem they should be able to fly into the detector without 
being forced into the hydrodynamic evolution. 

Fig. 2] illustrates the kr dependence of the HBT radii 
under various freeze-out conditions, which may be di- 
vided into two parts: 1) without HR and 2) with HR 
after the hydrodynamic phase. Here, without HR (lines) 
means that the evolution is stopped immediately after 
the Cooper-Frye freeze-out from the hydrodynamic phase 
(GF and 5 ■ £o are adopted as default hydro-freeze-out 
criteria), with instantaneous resonance decays. The ob- 
served size of the pion source is small at this hydrody- 
namic freeze-out. In previous investigations [2l| it has 
been found that binary baryon-meson and meson-meson 
collisions still frequently happen after the hydrodynamic 
freeze-out. In baryon-meson reactions, the most abun- 
dant interactions are the excitation and the decay of the 
A resonance (i.e. A ttN), while in meson- meson colli- 
sions, the p ^ 7T7T process is dominant. A large number 
of these final hadron interactions in which pions are in- 
volved contribute significantly to the final HBT radii of 
pions in this model. 

Let us therefore explore if the finally observed HBT 
radii do depend on the transition criterion from hydro- 
dynamics to the transport model. The full hybrid model 
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FIG. 4: kx dependence of the HBT radii Rl, Ro, and Rs (at 
midrapidity) for central HICs at SPS energies (-Eiab = 20A, 
30 A, 4(L4, 80Aand 158,4 GeV). The NA49 data are indicated 
by solid stars [48| . The HG-EoS is employed in all calculations 
but under different freeze-out conditions: 1) without HR, cal- 
culations with default hydro- freeze-out criteria (GF, 5 -so) are 
depicted by lines. 2) full hybrid model calculations with two 
different cuts of the energy density (default and 4 • eo) for the 
GF and with the default cut of the energy density for the IF 
are depicted by dashed lines with symbols. 



calculations (dashed lines with open symbols) are shown 
with two different cuts of the energy density (5 ■ £o as 
default and 4 ■ e ) for the GF and with the energy den- 
sity cut 5 • eo for the IF. It is found that the final state 
hadronic interactions are sufficient that the effects of dif- 
ferent treatments of the hydrodynamic freeze-out on the 
final HBT radii are almost totally washed out in all di- 
rections and at all investigated energies. 

Let us finally explore the dependence of the Ro/Rs 
ratio on the different EoS and freeze-out prescriptions. 
This ratio was expected to be sensitive to the duration 
time of the homogeneity region. In Fig. [5] the excita- 
tion function of the Ro / Rs ratio with the different EoS 
(lines with solid symbols) and freeze-out prescriptions 
(dashed lines with open symbols) are shown. The kx bin 
200 — 300MeV/c is chosen. The result for the pure cas- 
cade calculation is also shown as a baseline (dotted line). 
It is seen clearly that the Ro / Rs ratio is sensitive to the 
EoS, but not to the various hydrodynamic freeze-out pre- 
scriptions when including HR (shown as open triangles 
and open inverted triangles) as it has already been im- 
plied from the results of the HBT radii shown in Figs. [5] 
and |U With increasing latent heat which corresponds to 
the softness of EoS implied from Fig.Q] the Ro/Rs ratio 
is increased. The "excessively" large latent heat in BM- 
EoS results in a long duration time of the pion source and 
hence a large Ro/Rs ratio. Although the overall height 
is largely overpredicted by the BM-EoS, the qualitative 
behaviour of the data (with a maximal lifetime at beam 
energies around 40 — 100/1 GeV) is well reproduced. In 
addition, the "peak structure" is less pronounced than 
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FIG. 5: Excitation function of the Ro/Rs ratio at fcy = 
250MeV/c. The NA49 data are indicated by solid stars [i^ . 
UrQMD cascade calculation is shown by dotted line. Hybrid 
model calculations with EoS of HG, BM, and CH and with 
HR are shown by lines with solid symbols (the default hydro- 
freeze-out criteria (GF, 5-£o) are used). In the HG-EoS mode, 
various criteria of freeze-out, (IF, 5 • £0) and (GF, 4 ■ £0) with 
HR, and the default (GF, 5 • £0) without HR are shown by 
dashed lines with open symbols. 



in previous predictions [111 ], due to the different initial 
state and seems to provide a more reasonable estimate 
of the magnitude of the lifetime enhancement. The chi- 
ral EoS CH exhibits a lower Ro/Rs ratio because the 
first-order phase transition is less pronounced. The cal- 
culation with HG mode (line with solid squares) leads to 
the smallest Ro/Rs ratio due to the most stiffest EoS 
among the three ones. The result of the cascade calcu- 
lation lies in between the CH and the BM modes, which 
implies a relatively soft EoS. It can be understood since 
in the pure UrQMD model the new particle production 
is treated cither as a resonance decay or a fragmentation 
of the string, which introduces a finite lifetime and hence 
leads to a softer EoS. After considering the mean field 
potentials for both confined and "pre-formed" particles 
H3, > which gives a strong repulsion at the early stage, 
the Ro/Rs ratio was seen to decrease in line with results 
obtained here. 

For the full hybrid model calculation the different 
freeze-out prescriptions do not affect the final results 
(when comparing the results by dashed lines with open 
triangles with that by the line with solid squares). The 
calculation with HG-EoS but without HR (dashed line 
with open diamonds) seems to provide better description 



of the data, but, seems clearly unphysical to the authors 
as a solution to the duration time problem. 



V. SUMMARY AND OUTLOOK 

In summary, the HBT correlation of the negatively 
charged pion source created in central Pb+Pb collisions 
at SPS energies was investigated with a hybrid model 
that incorporates a (3+l)d hydrodynamic evolution in 
the UrQMD transport approach. We explored different 
settings, one where the EoS was varied without changing 
the initial conditions and the freeze-out prescription and 
another where the EoS was fixed and the treatment of 
the freeze-out was changed. We presented a systematic 
investigation of these effects on the HBT radii. It was 
found that the latent heat influences the emission of par- 
ticles visibly and hence the HBT radii of the pion source. 
While the final rescatterings result in an independence of 
the calculated HBT-parameters from the transition cri- 
terion, they do not improve the quantitative agreement 
with the experimental data. The details of the hydro- 
freeze-out prescription do not affect the HBT radii as 
well as the Ro / Rs ratio as long as the HR in the subse- 
quent hadronic transport model were taken into account. 
Overall, the HBT data seem to favor a stiff EoS [5l| , but 
one should also keep in mind that viscosity effects are 
neglected during the hydrodynamic stage and that the 
particle emission from the early stages should be han- 
dled more carefully. 

In the future, bulk and shear viscosity will be further 
considered for the hydrodynamic phase, and the non- 
equilibrated particle emission should be treated more pre- 
cisely. 
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